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Abstract. We give a rigorous derivation of the general-relativistic formula for the two-way 
Doppler tracking of a spacecraft in Friedmann-Lemaitre-Robertson- Walker (FLRW) and in 
McVittie spacetimes. The leading order corrections of the so-determined acceleration to the 
Newtonian acceleration are due to special-relativistic effects and cosmological expansion. The 
latter, although linear in the Hubble constant, is negligible in typical applications within the 
Solar System. 

PACS numbers: 95.30.Sf, 98.80.-k, 04.20.Cv 

1. Introduction 

Doppler measurements are very well studied in the case of weak-field approximation of the 
gravitational field of (arbitrary moving) isolated sources [1]. However, a similarly careful 
analysis of the cosmological effects on Doppler measurements is, in our opinion, still lacking. 
Considerations in this direction are contained in [2| for FLRW universes and in [3] for the 
Schwarzschild-de Sitter case. The aim of our paper is to give a rigorous derivation of the 
two-way Doppler formula in relevant cosmological spacetimes (FLRW and McVittie). 

The analysis of the radio Doppler tracking data from the Pioneer 10 and Pioneer 11 
spacecrafts yield an anomalous inward pointing acceleration of magnitude 8.5 x 10 -10 m s~ 2 ; 
see [4| and [5|. This magnitude is comparable to He = 7 x 10~ 10 ms~ 2 , where we set 
H = 72 kms" 1 Mpc -1 (see [6| for recent figures). This somewhat surprising coincidence 
invited speculations as to a possible cosmological origin of the 'Pioneer Anomaly'. Whereas 
there seems to be no disagreement over the absence of a genuine dynamical influence of 
cosmological expansion on Solar System dynamics, opinions on possible kinematical effects 
are less unanimous. Some even seem to claim that the Pioneer Anomaly can be fully 
accounted for by such effects [7| (for a critical discussion, see 0). Note that all these 
speculations rest on the assumption that the cosmological expansion extends into the small 
region occupied by our Solar System, i.e. that the expansion is not screened by local 
inhomogeneities of mass abundance, like that given by our Galaxy and further our Solar 
System. Our results here imply that even in the unlikely case that such a screening does not 
take place, there is no effect due to cosmological expansion of the required order of magnitude. 
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2. Basic kinematical definitions 

In this section we recall some basic kinematical notions that can be found, for example, 
in (9). Let (M, g, V) represent spacetime, where M. is a four-dimensionaQ manifold (space- 
and time-orientable) with Lorentzian metric g and Levi-Civita connection V. We adopt 
the 'mostly minus' signature convention where the restriction of g to spacelike directions 
is negative definite. We use units in which c = 1. For a vector v in the tangent space T p M. 
of M. at p 6 M. we define \\v\\ p := y/\g(v, v) p \. A normalized vector is one for which 

IMI = i- ' 

An observer (worldline) is represented by a timelike (smooth) curve in M. which 
w.l.o.g. we think of as being parameterized with respect to arc length and future-directed. 
By an observer at p S M. we understand a timelike future-pointing normalized vector in 
T p A4. An observer field is a timelike future -pointing normalized vector field on some open 
subset of M. . Let p be a point in Ai and u an observer at p. We define two projectors on T p M. 
by Q u '■= u <B> u b , and P u := id — Q u , where u b := •) is the one-form corresponding 
to u via g. We recall that the projection properties are as follows: P and Q are self-adjoint 
w.r.t. g, Q 2 = Q, P 2 = P, PQ = QP = 0, and P + Q = id. The local rest space of u 
at p is then defined as lZ u (p) := P u (T p A4). It locally represents the directions of the events 
that are Einstein-synchronized with the observer u at p. 

The above projectors can be naturally extended to general tensor fields in the standard 
way (e.g. for 1-forms one defines Pu> := u> o P), which we will denote by the same symbols. 
A tensor field T is called spatial w.r.t. u (at p) iff P U T = T (at p). Given an observer u at p 
we split g into two degenerate, positive semi-definite metrics 

t u := Q u g , (1) 

(T U ■= P U g , (2) 

which measure eigentime- and eigenspace-intervals with respect to u. 

Let now u and v be two observers at p. The relative velocity ofv w.r.t. u is given by 

P-U.V 



P u {v)p 
and its modulus by 



(3) 



Mv) p := \\f3 u ( V )\\ p = 



y/r u (v,v) 



V 1 g(u,v)i ■ l4i 



Given a unit vector e in the rest space of u (i.e. P u e = e and ||e|| = 1), the relative velocity 
of v w.r.t. u at p in direction e is given by 



0SMp : = °-™(X»> e ) p = -9(X»,e) p 



-fl(e,«) 



g(u,v) 



(5) 



Note that the modulus of the relative velocity of w.r.t. u equals that of u w.r.t. v. (Recall 
that the former velocity is measured with clocks and rods at rest w.r.t. u and is represented 
by a vector in lZ u (p), whereas the latter is measured with clocks and rods at rest w.r.t. v 
and is represented by a vector in lZ v (p).) Clearly < (3 u (v) < 1 and —1 < j3 u (v) < 1. 
Moreover ||Q„f|| = g(u,v) = l/\/l — f3 u (v) 2 is just the 'gamma-factor' familiar from 
Special Relativity. 



| Actually, all formal considerations that follow are independent of the dimension of spacetime. 
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The last notion we need to introduce in this section is that of relative acceleration of a 
worldline w.r.t. an observer field. For this purpose, given a worldline 7 and an observer field 
u along it, we define 

:= HQ^tH" 1 ^ V^oP„ . (6) 

V" is defined on tensor fields along 7 and P u stands for the extension of the spatial projector 
on the respective tensor space (i.e. every free index is to be projected). On functions we 
have V" = ||Q M 7||~ 1 7. Definition <JSJ resembles that of a rescaled (with the inverse of the 
'gamma-factor' ||Q U 7||) Fermi-derivative along 7. Notice, however, that the projection is 
taken w.r.t. u and not w.r.t. 7. It is easy to see that V™ is a derivation on spatial (w.r.t. u) 
tensor fields along 7. Moreover, it can be checked that V" is compatible with the spatial 
metric, i.e. the metricity property V" cr u = holds. 

We can now define the relative acceleration of a worldline 7 w.r.t. an observer field u 
along 7 as follows: 

-M - *?/>.«) - ^^(jgL^iMr) . ,7, 

which is a spatial (w.r.t. u) vector field along 7. Clearly, the relative acceleration of a worldline 
7 w.r.t. itself vanishes identically, i.e. 0^(7) = 0, as it should. For later use we also note that 

V^/3 M ( 7 ) 2 = 2<7 u (/3 u (7),a„(7)) , (8) 

and 

V?PZ(i) = * u (<x u ( 7 ),e)+a u (l3 u (>y),V%e) , (9) 

for some spatial (w.r.t. u) unit vector-field e along 7. Here we just used the metricity property 
of V? and /3 U ( 7 ) 2 = <t„(/3 u ( 7 ), /3J 7 )). 



3. Light rays 

Consider the propagation of a monochromatic electromagnetic wave. In the geometric-optics 
approximation (i.e. for wave-lengths negligibly small w.r.t. a typical radius of curvature of 
the spacetime and w.r.t. a typical length over which amplitude, polarization, and frequency 
vary) it propagates along a lightlike geodesic, along which the wave-vector field, k, obeys 
g{k,k) = Oand V k k = 0. 

An observer u at p who receives a light signal with wave vector k measures the frequency 

tj p u \k)=g(u,k) p . (10) 

Likewise, a second observer v at p will measure the frequency u> p v \k) = g(v,k) p . The 
relation between the two measurements is easily computed using (0]i and (0: 



^( fc ) = 1 (1 - bHv)) 



(id 



where A; := P u k/\\P u k\\ is the normalized vector in the rest space of u at p that points in the 
direction of the light propagation. This is nothing but the general Doppler formula. 

Let the observer v at p carry a mirror (meaning that the mirror is at rest w.r.t. v) that is 
used to reflect back the light ray. This means that the frequency, as measured by v, does not 
change whereas the spatial projection of the wave vector k changes sign. In short, the process 
of reflection is given by 

.(after) I _ ft . (before) | _ p . (before) | no) 



On Doppler tracking in cosmological spacetimes 



4 



Another observer, say u, at p will see a frequency shift according to 

^(feCaftcr)) _ 2 !-/%(„)„ i 
w («)( fe (bcforc)) l-jSuWI 

as one can compute from (V2\ , (fTTT i. and (0). Note that here and henceforth the wave vector k 
that defines k at the reflection point p refers to k just before reflection. 

4. Two-way Doppler tracking in a homogeneous and isotropic universe 

We consider a FLRW cosmological model, given by the metric 

g = dt 2 - a 2 (t)h , (14) 
together with the geodesic cosmological observer-field 

u = d/dt. (15) 

In ( fT4l /i is a homogeneous and isotropic metric (hence of constant curvature) on the slices of 
constant cosmological time t. 

Let k be the wave-vector field along a light ray (worldline). The frequency, as measured 
by the cosmological observers u, will change along the ray according to the well-known 
cosmological red-shift relation (see e.g. [ 10 1) 

uj^ u \k) (aot) = constant along the light ray . (16) 

The two-way Doppler tracking (see figure [TJ now consists in the exchange of light (radar) 
signals between us, the observers on Earth who hypothetic ally moves along the cosmological 
flow of u, and another observer 7 (spacecraft). Schematically the tracking involves the 
following five processes: (i) emission of the signal at po, (ii) propagation from pg to 
pi, (iii) reflection at p\, (iv) propagation from p\ to p%, and, finally, (v) reception at p 2 - 
Accordingly, using the cosmological red-shift relation (fTBl between po and px, the reflection 
shift (Qj} at px, and again the cosmological red-shift relation ( [Tol l between px and P2, one gets 

v£\k) _ a(t( Po )) ( l-p£(v) Pl \ 

4:\k) a(t( P2 )) { i-p u (v)i ) ■ 

This is the formula for the two-way Doppler tracking in a FLRW spacetime. It relates 
the spacecraft's velocity (relative to u) with the observable quantities (emitted and received 
frequencies, emission and reception times) if the scale function a is known. To get a better 
feeling of the above relation we note that to linear order in the two quantities HAfaa and (3 it 
reads 

« 1 - 2fc(v) Pl - H At 20 , (18) 

where At 2 a := t{p2) ~ t(po) is the coordinate-time interval between the emission and the 
reception events, which, since we are moving along d/dt, equals the eigentime interval 
measured by us between the two events. H := a/a denotes the Hubble parameter, which 
in ( TT8l can be evaluated at any of the times in the interval [to, fa]. 

Now, while tracking a spacecraft by continuously emitting a radar signal with constant 
frequency w the whole construction is pushed forward in time. One should think of the 
worldlines of the spacecraft (7) and ours as given. The latter is taken to be an integral curve 
of the cosmological observer- field u = d/dt and the former is given by some equation of 
motion (plus initial conditions), which we do not need to specify here. Thus, the three points 
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us 



spacecraft 




Figure 1. Sketch of the two-way Doppler tracking: The observer (us) moves on an integral 
curve of the cosmological vector field u. At po h e sends an electromagnetic signal to the 
spacecraft which moves along the worldline 7. The signal is then reflected back from the 
spacecraft at pi and, finally, received at p2 by the observer, fc := .P„fc/||.P, 1 fc|| is the 
normalized spatial (w.r.t. it) wave vector of the infalling light ray. 



(events) po, pi, and P2 are uniquely determined by any one of them. The same holds for the 
respective times i; := t(pi). We can then choose to express the events pi as functions of the 
reception time t2- Doing this, the reflection time t\ and the emission time to become functions 
of ti as well. 

One of the major tasks for the two-way Doppler tracking is to determine the spacecraft's 
(spatial) acceleration. Since ( ITTb relates the frequency shift with the velocity, differentiation 
of ( fTTT i w.r.t. the reception time gives a relation between the frequency-shift rate on one hand, 
and the acceleration and velocity on the other. In the differentiation w.r.t. t^, one has to take 
care of the different time dependences of the quantities in question. For example, the received 
frequency is to be thought of as function of the reception time <2 and the spacecraft's velocity 
as a function of the reflection time t\, which, in turn, is to be thought of as function of t^. In 
order to differentiate (TTTb w.r.t. <2 we thus need to know the dependences of the emission and 
reflection times, to and t\, on the reception time t^. These dependences can be obtained from 
the fact that the pairs of events (po,px) and (pi,P2) are lightlike separated. For simplicity, 
we will specialize to the spatially flat case, which is compatible with current observational 
data [6 |. Introducing spherical coordinates, (r,d,ip), on the slices of constant t, the spatial 
metric reads h = dr 2 + r 2 (d6 2 + sin 2 9 dtp 2 ). Without loss of generality we may assume our 
worldline to be given by r = 0. The spacecraft's worldline is described by some functions 
(ti,ri, 0±, (fx) of the reflection time tx- Hence one obtains 




(19) 
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where the minus sign is because the light ray is 'inward' pointing. Differentiating ( fT9T > w.r.t. t 2 
and noting that /3 k (v) Pl — a(ii)ri/ii = a(ti)dn/dti (where fc = \\d/dr\\~ 1 d/dr and the 
overdot denotes differentiation w.r.t. the eigentime of 7), one gets 

dti a(ti) 
a(t 2 ) 



(20) 



ri(ti(t 2 )) 1 

dr 

c 



dr 



ri(*l(*2)) 



Similarly, differentiating 

Jto(t 2 ) 

w.r.t. t 2 and using (|20T > one gets 

dtp = a{t ) ( 1- f3 k (v) Pl 
dt 2 a(t 2 ) \l + /3*(v) Pl 

In passing we remark that d22b can be used to compute the variation of the light round- 
trip-time At 2 o as measured by the receiver (At 2 Q is just twice the radar distance between 
observer and spacecraft): 



(21) 



(22) 



dh 



-At 



20 



1 



a(t 2 



HAt 



20 



(23) 



.1 + /3S(»)pi. 

which clearly displays the expected contributions due to the spacecraft's motion and the 
cosmological expansion respectively. 

Returning to the derivation of the frequency-shift rate, we notice that d/dt\ = 
(d£i/dri) -1 d/dTi = V" on functions along the spacecraft's worldline 7, where t\ is its 
arc-length. Taking into account ([8]), (O, ( |20l ), and ( 1221 ). differentiation of ( fTTl ) w.r.t. t 2 gives 

1 dW2 

(jJq dt 2 



(24) 



Here, we suppressed the arguments and the indices u for the sake of readability and put 




L0 2 



CO. 



(u) 

P2 



(fc). This formula gives the exact relation between the observable frequency-shift 
rate (measured 'here') and the local kinematical properties of the spacecraft (defined 'there') 
in a spatially flat FLRW spacetime — provided the scale function a is known. 

In the special, but relevant, case where the spacecraft's motion is in direction of the line- 
of-sight ('radial' for short), meaning that both (3 and a. are collinear to fc, we have (3 = (3 k fc, 
and ol = a k fc, where a k := <r(ot, fc). Moreover, it is straightforward to check that V"fc = 0. 
The two-way Doppler-tracking formula d24T > now simplifies to 



1 duj 2 
ojq dt 2 



a(t ) a(ti) 
a(t 2 ) a(t 2 ) 



2a k (l + [3 k )- 3 
I a(t 2 ) a(t ) 




(25) 
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This formula is still exact if one restricts to radially moving spacecrafts. In a linear 
approximation in HAt 2 o it becomes 

L^l M _ 2 f l-3If^ (a fe (l + /3 fe )- 3 +^ fe 1 • (26) 

In all cases of interest one is also interested in a slow-motion approximation involving an 
expansion in (3. Retaining only terms linear in j3 and HAt we simply get (dc^/d^V^o ~ 
— 2a h . This is exactly the result one obtains from Newtonian physics. Leading order 
corrections are obtained by going to quadratic order in (3 and also keeping mixed terms (3HAt. 
The result is 

_L ^fl ~ _ 2 (a* f l _ 3(3 k - 3HAt 2 o/2) + H(3 k ) . (27) 
ujq dt2 L V /J 

At this point it is crucial to understand how a relates to the acceleration that appears on 
the left-hand side of Newton's equation of motion. The latter is derived as the weak-field and 
slow-motion approximation of the spatial part of the geodesic equation. Recall that 'spatial 
part' refers to the choice of the observer field, from which the establishment of a Newtonian 
equation always depends. 

In addition, in order to write down a Newtonian equation, we need to use spatial 
coordinates of direct metrical significance. This is achieved by using a rescaled radial 
coordinate, r*(r, i) := a(t)r such that r* is now the proper geodesic distance to the origin 
(us) within each time slice of constant t. The time coordinate remains unchanged, :— t, so 
that 

d/dt* = d/dt - H(t)r d/dr . (28) 
The observer field is now the normalization of that field: 

u* = ||0/0i*|| _1 0/0t„ . (29) 

It also results from applying a (spacetime dependent) boost of parameter — {3 u (u*) = — -ffr* 
in the local planes spanned by u and e r := \\d / dr\\~ 1 d / dr. Note also that m» is hypersurface 
orthogonal, though clearly not orthogonal to the slices of constant t. 

The spatial projection of the geodesic equation with respect to up to quadratic order 
in H At and linear order in /3 Ut (7) turns out to be 

This is the form of the equation of motion in which dynamical considerations are most 
conveniently addressed (see [8| and references therein). Hence in the above formulae for 
the Doppler tracking we need to express the relative accelerations and velocities defined 
w.r.t. u in terms of relative accelerations and velocities w.r.t. it*. In the FLRW case these 
transformations are given in [8]. For a general spacetime they can be found e.g. in J9]. 
However, for our discussion it is sufficient to consider the approximated Doppler tracking 
formula (l27l and within its level of approximation it makes no difference whether we refer 
spatial acceleration and velocities w.r.t. u or w.r.t. m». Hence (|27] | gives the sought-after 
relation between the two-way Doppler measurement and the kinematical quantities that enter 
the Newtonian equation. 

Looking at (1271 we see, among other corrections, an additional acceleration term Hc(3 k 
(reintroducing the speed of light c). Hence the sometimes alleged He acceleration term |7) 
is actually suppressed by a factor j3, as already pointed out in [2|, though we cannot agree 
with the derivation given there. In particular we conclude that, even if the universe did expand 
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down to scales of the Solar System, there would be no effect of the same order of magnitude 
as the Pioneer Anomaly. Besides, notice that the above mentioned correction points in outer 
direction and hence opposite to the Pioneer anomalous acceleration. One expects that this 
conclusion remains true in the presence of an isolated local inhomogeneity, such as a single 
star, since the contributions coming from the 'gravitational' red-shift of the star cancel out 
in the two-way process. Insofar as such an inhomogeneity can be modeled by the McVittie 
metric this is explicitly shown in the next section. 



5. Two-way Doppler tracking in a McVittie spacetime 

The so-called 'fiat' McVittie model with central mass m and scale-factor a(t) is given by the 
spacetime metric 



1 — TOo/2a(t)r 



2 / N 4 



d f _ i + a\t)(dr z +rMtf) (31) 



l + m /2a(t)r J \ 2a(t)r 

together with a cosmological observer-field 

u = \\d/dt\\- 1 d/dt (32) 

along which the cosmological matter, given by an ideal fluid with pressure, moves. 

The red-shift formula for a radial light ray in a McVittie spacetime is easily computed up 
to linear order in HAtio and too / ar to be 

(\\d/dt\\u>M) t 



\d/dt\\uj( u )) po 



1 - HAt 10 , (33) 



where Aiio := £(pi) — t{po) an d H := a/a as previously. In case H = this reduces to 
the familiar gravitational red-shift relation in Schwarzschild spacetime, as it should. By ( 1331 
one sees that the main formula (Tl\ for the two-way Doppler frequency-shift we derived in 
the FLRW case still remains valid for the McVittie case in linear order in H At and to / ar. 
The reason is simply that the gravitational contribution to the red-shift (coming from the 
factors || <3/c?t || sa 1 — mo/ar in ( 1331 ) vanishes in the two-way process within the considered 
approximation. 

As in the previous section, we wish to compute the Doppler frequency-shift rate 
measured by an observer (us) along u who exchanges light (radar) signals with a spacecraft 
moving along an arbitrary worldline 7 (see figure [T]). For this we need the generally valid 
relations ([H) and (O and the relations < f20b and d22i > which, fortunately, hold exactly also in 
the McVittie case, provided that central mass, observer, and spacecraft are aligned. The latter 
will be assumed henceforth. 

A difference to the FLRW case is that here t is no longer the eigentime along the 
observer's worldline. Rather, to get the measured frequency-shift rate, we have to differentiate 
( fTTI ) with u = d/dr2 = ||<9/9i||~ 1 d/di2- When differentiating the kinematical factor 
on the r.h.s. of (T% we proceed as before but now take into account that d/dti = 
(dii/dri) _1 d/dri = \\d/dt\\ V™ on functions along the spacecraft's worldline 7, where 
n is its arc-length. This leads to the formula 

<34) 

which is valid in linear order in HAt2o and too /ar and quadratic order in ft, and differs from 
( |27] | merely in the term containing mo (note that in (|27| | Ai2o = At2o). Formula (l34l shows 
that in the McVittie case, too, the acceleration correction term He is suppressed by a factor ft. 
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As for the FLRW case, we may express the relative acceleration and velocity defined 
w.r.t. u in terms of the relative acceleration and velocity w.r.t. it*. As previously, it* is the 
normalized future-directed vector field whose integral lines have constant proper radius r* 
(i.e. radial geodesic distanc^ on the slices of constant t) and constant angular coordinates. 

Up to now 7 was unspecified. Let us now assume 7 to be a geodesic. In leading order in 
H At and mo /r* and linear order in Ut (7) the geodesic equation then leads to the following 
expression for the spatial acceleration w.r.t. u*: 



Here we made use of the identity a/a = —qH 2 , where q is the deceleration parameter, 
which shows that a/ a is of quadratic order in H, given that q is not a large number. Formula 
d35l l is the form of the equation of motion in the next-to-Newtonian approximation in which 
dynamical considerations are most conveniently addressed (see again [8| and references 
therein). 

In the special case of purely radial motion, insertion of $35[ into (f34t gives a formula 
predicting the two-way Doppler-shift rate in linear order in HAt2q and mo /r* and quadratic 
order in fi u ,{i)'- 



Hence there are two corrections to the Newtonian contribution. One is proportional to H 
and stems from the cosmological expansion, the other is independent of H and of purely 
special-relativistic origin, as one might check explicitly. Their ratio is (up to a factor ^/S) 
given by the square of the ratio of r* to the geometric mean of the Schwarzschild radius mo 
and the Hubble radius 1/H. The latter is of the order of 10 23 km, so that its geometric mean 
with a Schwarzschild radius of 1 km is approximately given by 2400 AU. The ratio of the 
effects is therefore of the order 10~ 7 . Hence the cosmological contribution is negligible for 
any application in the Solar System as compared to the special-relativistic correction. For the 
Pioneer spacecrafts 10 and 11 we have a radial velocity of about 12 km s -1 . This amounts 
to a special relativistic correction of magnitude 10~ 4 times the Newtonian gravitational 
acceleration, in an outwardly pointing direction. This is of the same order of magnitude 
as the Pioneers' anomalous acceleration J4], but pointing in the opposite direction. 

6. Conclusion 

It is often heard that effects of cosmological expansion are at most of order H 2 . We have 
seen in fl35l > that this is true as far as deviations from the Newtonian equation of motion 
are concerned. However, there is a term linear in H that enters the formula for the two-way 
Doppler tracking, as seen in (|36*T >. even though the correction it gives is negligible as compared 
with the special relativistic one, due to its suppression with a factor of (3. 

Finally we point out that for direct comparison with actual measurements there are other 
corrections to our Doppler formulae that can easily be taken care of. For example, one needs 
to incorporate the motion of the Earth relative to the cosmological substrate with an additional 
Doppler factor (fTTT l at the emission and reception points po and p2- Also, one needs to take 
into account that in the actual tracking procedure signals are sent back with a fixed frequency 
translation factor (see e.g. formula (1) in 0). 

§ In 1 8 1, for the McVittie case, we take the areal radius as radial coordinate for the Newtonian equation just for 
computational convenience. In the considered approximation, and in the region of validity of the Newtonian equation, 
the areal radius equals the geodesic radial distance. 




(35) 




(36) 
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